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Abstract 

We prove a global in time existence theorem, for the initial value prob- 
I lem for the Einstein-Boltzmann system, with arbitrarily large initial data, 

■ in the homogeneous case, in a Bianchi Type I space-time. 

a^ 
o , 

lO . 1 Introduction 

O 

Q [ The Einstein equations are the basic equations of the General Relativity, the 

cry funded in 1916 by Albert Einstein and according to which: "the geometric 
}L{ '. structure of the space-time is determined by its material and energetic con- 

W)' tents". These equations link, on one hand, the geometry of the space-time, 

which is determined by a fundamental tool: the metric tensor, standing for the 
gravitational field, and, on the other hand, all the material and energetic con- 
^ . , tents, summarized by the stress-energy-matter tensor (we will call it the matter 

■ tensor), acting as the sources for the gravitational field. Solving the Einstein 
equations is determining both the gravitational field and its sources. Naturally, 
the problem varies with the matter contents, and this is why, in the mathe- 
matical study of General Relativity, one of the main problems is to solve the 
Einstein equations coupled to various field equations, in order to determine both 
the gravitational field and its sources in different cases. A fundamental problem 
today is to establish the existence and to give the properties of global solutions 
to such coupled systems. We are interested in global dynamics of relativistic 
kinetic matter, and our goal is to prove such results in that domain. In the case 
of collisionless matter, the governing system is the Einstein- Vlasov system, in 
the pure gravitational case, or this system coupled to other field equations, if 
these fields are involved in the sources of the Einstein equations. In the colli- 
sionless case, several authors proved global results, see ^Tj, |2J for reviews, 
|12|. [22 and ^Uj for scalar matter fields, also see ^^li for the Einstein- 
Vlasov system with cosmological constant, which turns out to be a useful tool 
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for the proof of the fact that, the expansion of the universe is accelerating, see 
[TB] for more details on this question. 

Now, in the case of collisional matter, the Einstein- Vlasov system is replaced 
by the Einstein-Boltzmann system, which seems to be the best approximation 
available and that describes the case of instantaneous, binary and elastic col- 
lisions. In contrast with the abundance of papers in the coUisionless case, the 
literature seems very poor in the collisional case. If, due to its importance 
in collisional kinetic theory, several authors proved global results for the single 
Boltzmann equation, see |H1j 0) [3 for the non-relativistic case, and Pj for the 
full relativistic case, very few authors studied the Einstein-Boltzmann system, 
see ^ for a local existence theorem. It then seems interesting to us, to extend to 
the collisional case some global results obtained in the coUisionless case. This 
was certainly the objective of the author in d]i in which he studied the 
existence of global solutions for the Einstein-Boltzmann system. Unfortunately, 
several points of the work are far from clear, such as, the use of a formulation 
which is valid only for the non-relativistic Boltzmann equation, or, concerning 
the Einstein equations, to abandon the evolution equations, which are really 
relevant in the spatially homogeneous case he considers, to concentrate on the 
constraints equations, that reduce, as we will see, to a simple question of choice 
of the initial data. In [21] , the authors prove a recent result on the global exis- 
tence of solutions for the spatially homogeneous Einstein-Boltzmann system on 
a Robertson- Walker space-time, in the case of a strictly positive Cosmological 
constant A > 0. 

In this papers, we study the collisional evolution of a kind of uncharged 
massive particles, under the only influence of their common gravitational field, 
represented by the metric tensor we denote by 17, and which is a function of the 
position of the particles, and whose components gap, sometimes called 
"gravitational potentials", are subject to the Einstein equations. Here, the 
sources of the Einstein equations are generated by the only matter contents 
which are the massive particles, through the matter tensor which is a symmet- 
ric 2-tensor, we denote T^p- The particles are statistically described in terms 
of their distribution function, we denote by /, and which is a non-negative real- 
valued function of both the position and the momentum of particles. The scalar 
function / which is physically interpreted as the "probability of the presence 
density" of the massive particles, during their collisional evolution, is subject 
to the Boltzmann equation, defined by a non-linear operator called the " col- 
lision operator" . In the binary and elastic scheme due to Lichnerowicz and 
Chernikov (1940), we adopt, at a given position, only 2 particles collide each 
other, in an instantaneous shock affecting only the momentum of each particle, 
only the sum of the 2 momenta being preserved. We then study the coupled 
Einstein-Boltzmann system in (g, /). The system is coupled in the sense that 
each unknown appears in the equation determining the other unknown: / which 
is subject to the Boltzmann equation generates the sources of the Einstein equa- 
tions through the matter tensor Tap, whereas the metric tensor g which is sub- 
ject to the Einstein equations is in both sides of the Boltzmann equation in /. If 
the particles were charged, the above system would be coupled to the Maxwell 
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equations that would describe the electromagnetic effects; this is our project in 
a future work. 

We now specify the geometric framework, i.e. the kind of space-time we are 
looking for in the present work. An important part of the General Relativity 
is Cosmology, which is the study of the universe at large scale; in such a view- 
point, even cluster of galaxies are assimilable to "particles". A. Einstein and 
W. de Sitter introduced cosmological models in 1917. 

Let us point out that the Einstein equations are overdetermined and physi- 
cally meaning symmetriy assumptions reduce the number of unknown. The 
surface symmetry assumptions which are the spherical, plane and hyperbolic 
symmetries assumptions constitute the major part of the models studied in the 
works we quote in the references. Robertson and Walker showed in 1944 that 
" Exact spherical symmetry about every point would imply that the universe 
is spatially homogeneous" see [H], p- 135. In the Robertson- Walker models, 
the spatial geometry has constant curvature which is positive, zero and neg- 
ative respectively. We look for a spatially homogeneous, locally rotationally 
symmetric (LRS) Bianchi type I space-time, which is a direct generalization of 
the Robertson- Walker space-time with zero curvature and which is known, in 
Cosmology, to be the basic model for the study of the expanding universe. In 
the model we are looking for, the metric tensor g has only two unknown com- 
ponents a and b, and the spatial homogeneity means that a and b depend only 
on the time t, and the distribution function / depends only on the time t and 
the momentum p of the particles. The study of the Einstein-Boltzmann system 
then turns out to the determination of the triplet of scalar functions (a, 5, /), a 
and b being strictly positive functions. 

We now sketch the strategy we adopt to prove the global in time existence of a 
solution (a, b, f) of the initial value problem for the coupled Einstein-Boltzmann 
system, for arbitrarily large initial data (00,601/0) at the initial time t = 0. In 
the spatially homogeneous case we consider, the Einstein equations are a system 
of three non-linear o.d.e in a and b. The Boltzmann equation is a non-linear 
first order p.d.e for the distribution function /. 

In the first step, we suppose a and b given with the only assumption to be 
bounded away from zero; and we give, following Glassey R.T. in the correct 
formulation of the relativistic Boltzmann equation in / on a Bianchi type I space- 
time. We then prove that, on any bounded time interval I — [to,to -\- T], with 
to e M+, T e R*^^, the initial value problem for the Boltzmann equation, with 
initial data ftg € LKR."^), ft„ > a.e; has a unique solution / G C[I; L\iM.^)\, 
f{t) > a.e, yt € /, where Ll{M.^) is a weighted subspace of L^(R'^), whose 
weight , is imposed by the manner in which the sources terms Tq,^ of the Ein- 
stein equations depend on /. We follow the method developed in 23 to prove 
the global existence of the solution / G C[0, -l-oo; L^(R'^)] for the initial value 
problem for the Boltzmann equation, but here the weighted norm of Ll{M.^) 
could allow us to prove the existence theorem only on bounded time intervals 
I — [tQ,tQ + T], but this was enough for the coupling with the Einstein equations 
and to obtain the global existence by another method. 

In a second step, we suppose / given in C[I; L2{M.'^)] and we consider the Ein- 



3 



stein equations in a and 6, that split into the constraint equations and the evo- 
lution equations. The constraint equations contain the momentum constraints 
which are automatically satisfied in the homogeneous case we consider, and the 
Hamiltonian constraint which reduces to a question of choice for the initial data. 
The main problem is then to solve the evolution equations which are a system 
of two non-linear o.d.e in a and h. We set, following Rendall, A.D., and Uggla, 
C. in [m 

H=-— 1 y 

3 ' ^ a-2 + 26-2 + 1' 62 + 2a2' + trkb 

in which k is the second fundamental form on the space-like 3-manifolds that 
foliate the space-time in the 3+1 formulation of the Einstein equations. We 
then prove that the Einstein evolution system in a and b is equivalent to a first 
order system in (iJ, z, s, for which standard theorems apply. From there, 
we deduce the existence of the solutions a, h of the initial value problem for 
the Einstein evolution equations, on any interval / = [0,T[, T > 0, which are 
increasing function of t, and hence bounded from below, satisfying consequently 
the assumptions of the first step. 

In a third step, we use the change of variables indicated above together with 
the characteristic method for the first order p.d.e that reduces the Boltzmann 
equation to a non linear o.d.e in / with values in the Banach space L\{S?), 
to prove the local existence of solutions for the coupled Einstein-Boltzmann 
system by applying standard results. 

In a fourth step, we use the results of the two first steps to construct an 
appropriate functional framework, allowing us to prove, using the fixed point 
theorem that, the local solution of the coupled system, whose existence is es- 
tablished in the third step, is in fact, a global solution. 

Our method preserves the physical nature of the problem that imposes to the 
distribution function / to be a non-negative function, and this is why we had to 
consider spaces of functions / which are positive everywhere. Nowhere we had 
to require smallness assumptions on the initial data, which can, consequently, 
be taken arbitrarily large. 

An aspect of the results of the present paper is an extension of the results of 
j24| . to the case of the zero cosmological constant A = 0. Notice that the present 
paper follows papers [23 and [23] and the proofs for the Boltzmann equation 
are similar, but profound differences arise with |24j on the Einstein equations. 

This paper is organized as follows: 
In section 2, we introduce the Einstein equations, the Boltzmann equation and 
the coupled Einstein-Boltzmann system, on a Bianchi type I space-time. 
In section 3, we study the Boltzmann equation in /. 
In section 4, we study the Einstein equations in a and h. 

In section 5, we prove the existence of local solutions to the coupled Einstein- 
Boltzmann system. 

In section 6, we prove the existence of global solutions, to the coupled Einstein- 
Boltzmann system. 
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2 The Boltzmann equation, the Einstein equa- 
tions and the Einstein-Boltzmann system on 
a LRS Bianchi type I space-time. 

2.1 Notations and function spaces 

A greek index varies from to 3 and a Latin index from 1 to 3, unless otherwise 
specified. We adopt the Einstein summation convention aab" = '^aab". We 

a 

consider a locally rotationally symmetric (LRS) Bianchi type I space-time de- 
noted (M^, g) where, for x = {x") = (x^jX*) G M^, x^ —t denotes the time and 
X = (a;*) the space; g stands for the metric tensor with signature (-, +, +) 
that can be written: 

g = -dt^ + a^{t){dx^f + h'^{t)[{dx^f + {dx^f] (2.1) 

in which a and h are strictly positive functions of t. We consider the collisional 
evolution of a kind of uncharged massive relativistic particles in the time 
oriented space-time (R^,g). The particles are statistically described by their 
distribution function we denote /, which is a non-negative real-valued func- 
tion of both the position (x") and the 4- momentum p = (p") of the particles, 
and that coordinatize the tangent bundle T(M*) i.c: 

/ : r(M'*) ~ M"* X M'' ^ IR+, (x",p") ^ /(x",p") G M+ 
For p = (p*), q ~ (g') g R"^, we set, as usual \p\ — define the 

i 

scalar product: 

<p.q>=a'^p\^ +b^[p^q'^ +p\^] (2.2) 

We suppose the rest mass to > of the particles normalized to unity, i.e to = 
1. The relativistic particles are then required to move on the future sheet of 
the mass-shell, whose equation is g{p,p) = —1; from which we deduce, using 
expression H2.1|l of g: 

P° = v/l + a2(pl)2+62[(p2)2 + (p3)2] (3.3) 

where the choice of p'' > symbolizes the fact that, in a natural manner, 
particles eject towards the future. I|2.3|l shows that, in fact, / is defined on the 
7-dimensional subbundle of T(K'') coordinatized by (a;"), (p^). In the spatially 
homogeneous case we consider, / depends only on t and p = (p*). 
The framework we will refer to is the subspace of L^(]R'^), denote L2{M.'^) and 
defined by : 

LliW") ^ {f e L\R^),\\f\\ -.^ f y^lTW\fip)\dp<+^} (2.4) 

is a norm on Ll{M.^) and (L2(M''), ||.||) is a Banach space. We set for r an 
arbitrary strictly positive real number; 

Xr = {feLl{R'), />0 a.e, \\f\\<r} (2.5) 
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One verifies easily that, endowed with the metric induced by the norm . , JY"^ 
is a complete and connected metric subspace of {Ll{M.'^), ||.||). Let I be an real 
interval; we set: 

C[I; Ll{R^)] ^{f : I LI{M.^), f continuous and bounded } 

endowed with the norm: 

Ill/Ill = ^HI/WII (2.6) 
tei 

C[I; Ll{M.^)] is a Banach space. Xr being defined by (|2.5|l : we set: 

C[I;Xr]={feC[I;Ll{R% f{t)eXr, Vt e 1} (2.7) 

Endowed with the metric induced by the norm |||.||| defined by H2.6f) . C[I;Xr] 
is a complete metric subspace of (^C[I; L2{R^)], |||-|||)- 

2.2 The Boltzmann equation 

In its general form, the Boltzmann equation in /, on the curved space-time 
{W^, g) can be written: 

^'"^-r»^^ = Q(/,/) (2.8) 

Where (r"^) are the Christoffel symbols of g, and Q is a non- linear integral 
operator called the " Collision Operator" . We specify this operator in details in 
next section. Now, since / depends only on t and (p*), (|2.8|l writes : 

p'%-Vp'§^ = Q{fJ) (2.9) 

We now express the Christoffel symbols in the case of the metric g defined by 
H2.1|l . Their general expression is: 

r^/^ = Ig^'^idc.g^^ + dpg^^ - d^g^^p] (2.10) 

in which {g^^) denotes the inverse matrix of {g\^j)\ We have F^^ = F^^. The 
definition (|2.1|) of g gives at once: 

„00 _ „ _ 1 . „ — „2. „ _ „ — ^2. „11 _ „-2. „22 _ „33 _ 1,-2 

9 — .900 — ~i, 5ii — a I .922 —533 — 0, 9 —a , g —g — o 
50* = 5" 0; 5y = g'^ =0 for i^j 

(2.11) 

A direct computation, using H2.10|l and (|2.11() . gives, with a = ^: 

{T^i — da; = F33 = bb; FJq = ^; Fjq = F3Q = |; Fqq = 0; 

FO^ = for a^(i- F^- = 

(2.12) 
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The Boltzmann equation (|2.1U|I then writes, using (|2.12ll : 

in which p*^ is given by (|2.3() . (|2.13|) is a non-hnear p.d.e in / we study in next 
section. 

2.3 The Einstein Equations 

The Einstein equations 'm g — {gap) can be written: 

i?a/3-^i?ga/3-87rGr,/3 (2.14) 

in which: 

Rafj is the Ricci tensor , contracted of the curvature tensor of g we specify 
beUow; 

R = g°'^Rap = Ra is the scalar curvature; 

Taf3 is the matter tensor that represents the matter contents, namely the massive 
particles in our case. Tap is generated by the distribution function / of the 
particles by: 

T-^t) = / tl^m^dp (2.15) 

where \g\ is the determinant of g; H2.1|l gives \g\^ ~ ab^. Recall that / is a 
function of t and p — (p*); hence T"'^ is a function of t. 
G is the universal gravitational constant. We take G — I. 

The contraction of the Bianchi identities gives the identities VqS'"'' = 0, 
where Sap = Rap — \R dap is the Einstein tensor. The Einstein equations H2.14|l 
then imply that the tensor T°'P must satisfy the four relations VqT"^ = called 
the conservation laws. But it is proved in 2 that these laws are satisfied for all 
solution / of the Boltzmann equation. We have: 



Rap - R a,xp 
where 

^li',ap = ^oX^p — dpT^^ + ^ta^'^P ~ ^ip^lla 



(2.16) 



in which F^^ is given by H2.10|l with the only non zero components given by 
H2.12|l . H2.12|l and H2.16|l show that the Einstein equations l|2.14|l are a system 
of second order non-linear o.d.e in a and b. In order to have things fresh in mind 
at the adequate moment, we leave the expression of (|2.14ll in term of a and b to 
section 4, which is devoted to the study of the Einstein equations. 

2.4 The coupled Einstein-Boltzmann system 

H2.14|l - H2.13l) is the Einstein-Boltzmann system we study. The system is coupled 
in the sense that, / which is subject to the Boltzmann equation H2.13|l generates 
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through H2.15|l . the sources terms Ta/s of the Einstein equations; it appears 
clearly that the coefficients of the Boltzmann equation 12.13|l in / depend on 
a and b which are subject to the Einstein equations (|2.14|l . and we will see, in 
next section, that a and b also appear in the collision operator Q, which is the 
r.h.s of the Boltzmann equation (|2.13|) . 

3 Existence Theorem for the Boltzmann Equa- 
tion 

In this section, we suppose that the components a and b of the metric tensor g are 
given and fixed, and we prove an existence theorem for the initial values problem 
for the Boltzmann equation H2.13|l . on every bounded interval I — [to,to + T] 
where G K.+ , T € M^. We begin by specifying the collision operator Q in 

3.1 The Collision Operator 

In the instantaneous, binary and elastic scheme, due to Lichnerowicz and Chernikov 
we consider, at a given position {t,x), only 2 particles collide instantaneously, 
without destroying each other, the collision affecting only the momentum of 
each particle that changes after the shock, only the sum of the 2 momenta 
being preserved, following the scheme: 




p + q = p' + q' 

The collision operator Q is defined, using functions /, g on M.^, p, q standing for 
the momenta before the collision and p' , q' the momenta after the collision by: 

Qif,9)^Q^if,9)-Q-if,9) (3.1) 

where: 

Q+{f,g){p)= f f{p')g{q')A{a,b,p,q,p',q')di, (3.2) 

Q-{f,gKp)= I / f{p)g{q)A{aAp,q,p'.q')d^ (3.3) 

whose elements we now introduce step by step, specifying properties and hy- 
potheses we adopt: 
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1) S'^ is the unit sphere of K^, whose volume element is denoted dw. 

2) A is a non-negative real-valued regular function of all its arguments, called 
the collision kernel or the cross-section of the collisions, on which 
we require the following boundedness, symmetry and Lipschitz continuity 
assumptions: 

< A{a, 6, p, q, p', q') < Co (3.4) 

A{a, b, p, q, p', q') = A{a, b, p', q', p, q) (3.5) 

A{a, b, p, q, p', q') = A(a, 6, q, p, q' , p') (3.6) 

\A{ai, b, p, q, p', q') ~ A(a2, b, p, q, p', q')\ < ko\ai - aa] (3.7) 

\A{a, bi, p, q, p', q') - A{a, 62, P, q, p', q')\ < ko\bi - 62I (3.8) 
where Cq and ko are strictly positive constants. 

3) The conservation law p + q = p' + q' splits into: 

/ + g"=P° + g" (3.9) 

p + q = p' + q' (3.10) 
and (|3.9|) shows, using H2.3|l . the conservation of the quantity: 



(3.11) 

called the elementary energy of the unit rest mass particles; we can inter- 
pret (|3.1U|) by setting, following Glassey, R. T., in |11| : 



\p> ^p + R{p, q,uj)uj 
\q' — q — R[p, q, iL!)uj ; uj £ S' 



(3.12) 



in which R{p,q^Lo) is a real-valued function. We prove, by a straightfor- 
ward computation (using H2.3|l to express p' , q' in terms oi p', q', and 
next (|3.12() to express p', p' in terms of p, p, lu, R, that equation (|3.9|l 
leads to a quadratic equation in R, that solves to give the only non trivial 
solution: 

e-^ - [< LO.{p + q)Y 

in which p — -tt, e is given by H3.11|) . and <, > is the scalar product defined 
by (|2.2|) . Another direct computation gives, using the properties of the 
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determinants, that the Jacobian of the change of variables {p, q) {p', q') 
in R3 X M^^ defined by !^J^ is: 

d{p,q) p°q° 

^A.1A\ shows, using once more (I2.3|l and the imphcit functions theorem, that the 
change of variables 13.12|l is invertible and also allows to compute p, q in terms 
of p' , q' . 

Finally, formulae and (|3.12|) show that the functions to integrate in (|3.2|l 
and 1)3. 3|l completely express in terms of p, w, and the integrations with respect 
to q and u leave functions Q~^{f, g), Q^{f, g) in terms of the single variable p. In 
practice, we will consider functions / : R x — > R, that induce, for t fixed in R, 
functions f{t) on R'^, defined by f{t){p) — f{t,p). Even in such cases, in order 
to simplify notations, we will explicit the dependence on t only if necessary. 

Remark 3.1 1) Formulae 1^3.1S\I and {3.14}) are generalizations to the case 
of the LRS Bianchi type I space-time, of analogous formulae established 
by Glassey, R. T, in tllf . in the case of the Minkowski space-time, which 
corresponds to the case a{t) = b{t) = 1 in the metric defined by 1^2. 1]) . in 
which case the scalar product <,> defined by \2.!l^) reduces to the usual 
scalar product in R^. 

2) In ujq = \g\^^ — ob^^ is a Leray form, which is in fact, the 
canonical volume element in the momenta space. 

3) A — fce^° ^'^ ^'p' "'"^'l , k > 0, is a simple example of collision 
kernel satisfying assumptions {3.4^ , 1^3.5}) . ^3.6]) . \3. ?| ) and ^3.8}) . 

3.2 Resolution of the Boltzmann equation 

We consider the Boltzmann equation (2.13) on / = [to, tg + T], to G R+, T > 0. 

The functions a and b are supposed to be defined on /. Equation l|2.13|l is 
a first order p.d.e in / and its resolution is equivalent to the resolution of the 
associated characteristic system, which can be written, taking t as parameter: 

dp^ a 1 dn^ b 9 dp^ b n df 1 , , , 

ir--'y-^ -ir-~'bP-^ 1--^/= i-^w,/) (3.15) 

We solve the initial value problem on / = [to, to + T], with initial data: 

p'(to)-y^i = 1,2,3; /(to) = /to (3.16) 
The equations in p = (p*) solve directly to give, setting y = (y*) £ M.^; 



^ , r./ N ^ / a^(^o) b'^ih) &^(to) 
p{to+t,y) = D{t)y, with D{t) = Diag / ° ^ u2u \v l^IT^ 

\ a-^ (to + t) r (to + t) r [to + tj 



, t e [0 
(3.17) 
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which shows that, the initial value problem for the Boltzmann equation (2.13) 
is finally equivalent to the following integral equation in which, for simplicity, 
we conserve the notation p, which stands this time for any independent variable 
in M3. 

f{to+t,p)^ftM+ f"^ ^QifJ)is,P)ds te[0,T] (3.18) 

Jto P 

We prove: 

Theorem 3.1 Let a and b be strictly positive continuous functions oft, such 
that a{t) > |, b{t) > |, whenever a and b are defined . Let ftg S L2{M.^), 
/to ^ 0; '^•fi; fl^*^ € M*^ such that r > \\fta\\- 

Then, the initial value problem for the Boltzmznn equation on [ig, to + T], with 
initial data fto, has a unique solution f S C[to, ^r] • Moreover, f satisfies 

the estimation: 

Sup 11/(011 < ll/to II (3.19) 

ie[to,to+T] 

Proof: Theorem 13. II is a direct consequence of the following result: 
Proposition 3.1 Assume hypotheses of theorem \3.1\ on: a, b, ftg andr. 

1) There exists an integer no{r) > 1 such that, for every integer n > no{r) 
and for every v S Xr, the equation: 

\/nu -Q{u,u) = v (3.20) 

has a unique solution Un G Xr- 

2) Let n eN, n> no{r) 

i) For every u g Xr , define R(n, Q)u to be the unique element of Xr 
such that: 

^/nR(n, Q)u - 4t<3 Q)u, R(n, Q)u] = u (3.21) 

ii) Define operator Qn on Xr by: 

Qn{u,u) — n^/nR(n,Q)u — nu (3.22) 

Then 

a) The integral equation 

f{to+t,p) = ftM + Qn{fJ){s,p)ds t e [0,T] (3.23) 

Jto 

has a unique solution /„ € C[tQ,tQ -'rT; Xr]- Moreover, f„ satis- 
fies the estimation: 

Sup ||/„(t)|l < ll/tJI (3.24) 

t£[to,to+T] 
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b) The sequence {fn) converges in C[to,tQ + T;Xr]to an element 
f € C[tQ,to + T; Xr], which is the unique solution of the integral 
equation |X7^. The solution f satisfies the estimation 1^3. iy\) . 

Proof of the proposition 13. ll 

As we pointed it out, the above result is similar to those of |^ and The 
proof of prop 13.11 above is analogous to the proof of prop. 3.1 of [23, and it 
follows the same lines as the proof of theorem 4.1 in We will emphasize 

only on points where differences arise with the present case, and, at the same 
time, we establish useful formulae for next steps of the present paper. 

A) Proof of point 1) of prop 3.1 
We use: 

Lemma 3.1 Let /, .g G L1(M3). then j,Q+ {f , g), ^Q- (f, g) e LUR^) and 



<C{t) 11/1111,9 11, 



<C{t)\\f WW g\\ (3.25) 



joQ+ifJ)-joQ+ig,g) 
1 



if, f) — 



1 



(5.5) 



<CW(ll/ll + l|g||)ll.f-5ll 

<CW(||/|| + ||.g||)||/-5ll 



(3.26) 



where 



^QifJ)-^Qi9,9) 



<c{t)iWfW + W .9 II) II /-5 II 



C{t) = Ca62(i) 
in which C > is an absolute constant. 



(3.27) 
(3.28) 



proof of lemma 3.1 

We deduce from l|3.9|l . using a > 1, > 1 and expression H2.3() of p° 



1 



V^TW < y 1 + ^ + < 2(/ + g°) = 2(p'0 + g'") 
Expression (|3.2|1 of Q^{f,g) then gives, since by H3.4|l . \A\ < Cq: 



-oQ+(/,5) 



- — — Q if'9) 



dp 



< 2ab^{t)Co 



\fip')Mq')\dpdqdiu 



We make the change of variable 13.1211 and we deduce from (|3.14|) that gives 

dpdq = dp'dq', and since by (jTSI) we have ^td^ = ^ + ^ < 2, that 

the above inequality yields: 



1 



?+(/,5) 



< ^TTab^{t)Co 



\f{p')\\9{q')Wdq' <Cab'Wf\\\\9\\ 
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follows the same way without 



with C = SttCq. The estimation of -^Q {f,g) 

change of variables and (|3.25|) follows. The inequalities H3.26|l are direct conse- 
quences of (|3.25|l and the bilinearity of (5+ and , which allows us to write, 
P standing for or ^Q": 

Pif, I) - P{9. g) = P{f, f-9) + Pif - 9, 9)- 

Finally, (|3.27|l is a consequence of (|3.2t)|) and Q ~ — . This completes 
the proof of the lemma 3.1H 

Now, since the function t ab'^{t) is positive and continuous on the line 
segment [to, to +T], there exists an absolute constant C{to,T) > such that 
C{t) < C{to,T), yt e [to, to +T], where C(t) is defined by ^X^ . We then 
obtain, from lemma 3.1, the same inequalities with an absolute constant, than 
the inequalities in proposition 3.1 in |23| . The proof of the point 1) of propo- 
sition 3.1 above is then exactly the same as the proof of proposition 3.2 in [23]. ■ 

B) Proof of the point 2) of prop 3.1 

We use, n(){r) being the integer defined in point 1) of proposition 3.1 above : 

Lemma 3.2 We have, for every integer n > no(r) and for every u G Xr 

\\V^R{n,Q)u\\^\u\\ (3.29) 

Proof of the lemma 3.2 

H3.29|l is a consequence of : 



/ QifJ)iP)dp = 0, yfeLi 

JR3 



(3.30) 



H3.3()|l is established exactly as formula (3.30) in the proof of lemma 3.5 in )24| . 
replacing everywhere in that proof, by ab^, and since we did, in the present 
paper, the same assumptions and we have the same properties as those used 
in pi]. Now deduce (|3.29() from H3.30|l exactly as in the proof of lemma 3.5 in 
|24|. choosing this time B — Diagia, b, b) instead of _B = Diagla, a, a) and using 
a > |: ^ > |: to conclude the proof of lemma 3.2 above. 

Now (|3.29|l is exactly the equality 13.10|l in proposition 3.3 in |^. We then 
prove exactly as for proposition 3.3 in 23 , that all the other relations hold in the 
present case. The proofs of points 2) a) and 2) b) of the above proposition 3.1 
are exactly the same as the proofs of proposition 4.1 and theorem 4.1 in |23| . 
just replacing [0, -|-oo[ by [to, to + T]. This completes the proof of proposition 3.1 
above, which gives directly theorem 3.1 ■ 



4 Existence Theorem of the Einstein Equations 

4.1 Expression and Reduction of the Einstein Equations 

We express the Einstein equations H2.14|l in terms of a and b. We have to 
compute the Ricci tensor Rai3 given by H2.16|l . The expression (2.12) of F^o 
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shows, using (2.16), that the only non-zero components of the Ricci tensor are 
the Raa and R22 = -Ras- Then, it will be enough to compute i?oo — R\ xq] 
Rii = R\ and R22 — R\ \2- "^^^ expression ()2.12fl of F^^ and formulae 
give: 



-^"0.20 — R Q. 



30 



-^"1,01 — R\^ii = 0; R\^21 — R\,31 — 0-0.- 



We then have: 



a 



b b •• „ 

RoQ = 2-; Rii = ad + 2aa-; R22 ^ R33 ^ bb + bb- + (by 

a b b a 



We can then compute the scalar curvature R to be: 

R = g"^Rci3 = 2 



a b ab 
- + 2- + 2- 
a b ab 



The Einstein equations (|2.14|) in a and b then take the reduced form: 

ab f bV , , 



-6^ 



a6 
a5 



SvrTn 



(4.2) 



(4.3) 



in which T^^ = gaxQiSfiT^^ is given in terms of / by (|2.15|) . 
In this paragraph, we suppose / fixed in C[/; Xr], where / = [0, T], T > 0, with 
/(O) = /o e i^K^). /o > a.e. and r > ||/o||, r fixed. We study the initial 
value problem for the non- linear second order system H4.1() - (I4.2() - (|4.3() . in a and 
6, with initial data oq, b^, do, &o; i-e 



a(0) = ao; 6(0) = 60; a(0) = do; 6(0) = 60 



(4.4) 
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4.2 Compatibility 

We notice that, if Sap = Rap — jdapR is the Einstein tensor, we have 

S'oi = 0; Sij =0 for j; S22 = S33 
then, the compatibiUty of the Einstein equations 12.14|l require that: 

To^^0; Ty ^O for i^j; T22 - T33 (4.5) 

But the matter tensor Tap = gaxgpfiT^^ is defined by H2.15|) in terms of the 
distribution function /. It then appears that (|4.5I) are in fact conditions to 
impose to /. we prove: 

Proposition 4.1 Let G be the sub-group of 0(3) whose elements are on the 
form: 

/ 6 \ 

M,(6i) = cos6l -sine ; = 1, 61 e M 
\ smB cos 6* / 

Assume the hypotheses of theorem 3.1 for t^ — Q. Assume in addition that /o is 
invariant by G, and that the collision kernel A satisfies: 

A{a, b, Mp, Mq, Mp', Mq') = A{a, b,p, q,p' , q'), Vp, q G M? , VAf e G (4.6) 

then 

1) The solution f of the integral equation ^y.lS\) on [0,T] satisfies: 

fit,Mp)^ fit,p) Vte[0,T], ypeR^ yMeG (4.7) 

2) The matter tensor Tap satisfies the conditions 

Proof of Proposition 4.1: Let AI € G 

1) One verifies easily, that the scalar product <, > defined by H2.2|l and, conse- 
quently, p^{p) defined by H2.3|l are invariant by G, which means 
< Mp,Mq >=< p^q p°{Mp) = p°{p), VM e G. Now (|XTH|l gives, since 
fo{Mp) = foip): 

fit,Mp)^foip)+ f ^Q{fJ)is,Mp)ds, Vte[0,T] (4.8) 
Jo P 

Next, definition (|3.1() . H3.2(l . H3.3|l of the collision operator Q gives: 

Q{fJ){s,Mp)^ / ^dq [f[s,p')f{s,q')~f{s,Mp)f{s,q)]A{aAMp,q,p\q')dLO 

(4.9) 

Let us set in H4.9|l q = Mqi; ut = Mloi. Then, formulae 13.12|l give, using 
expression H3.13|l of R and the invariance of the scalar product <,> by the 
subgroup G of 0(3), VM e G: 

r p' = Mp + R{Mp,q,uj)uj = Mp + R{Mp, Mqi, Muji)Mlji = M {p + R{p,qi,uJi)uji) 
\ q' ^ q - R{Mp,q,uj)uj = Mqi - R{Mp, Mqi, MoJi)MuJi = M (qi - R{p,qi,uJi)LUi) 
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So that p' = Mp\; q' = Mq\ where: 

p\ ^ p + R{p,qi,uJi)uJi; q\ ^ qi - R{p, qi,uji)uji. 

Then (|4.9|) impHes, using assumption H4.6I) on A, ~ qi, and the invariance of 
the volume elements dq, dto by G: 

Q{fJ){s, MP) = Q[f{s)oMJ{s)oM]{p) 

H4.8|l then writes: 

(/(t)oAf) (p) = /o(p) + / ^Q[fis)oM, f{s)oM]{p)ds 
Jo P 

which shows, by setting h{s) = f{s)oM, that ||/i(s)|| — (by setting 

q = Mp in the integral defining ||/(s)oAf ||) , and that, h and / are 2 solutions 
in C[Q,T; Xj.] of the integral equation (|3.18|) . The uniqueness theorem 3.1 then 
implies that h = f and the point 1 of proposition 4.1 is proved. 
2) Consider the expression H2.15|l of T"^ in which / satisfies (14.711 and observe 
that, if p = Mq where M G G, then = p°{Mq) = p°{q) q° 

(i) Set in H2.15|l a — 0, (3 — i with i = 1, 2, 3; now compute the integral using 
the change of variables p = 7\f_i(7r)g; where M^{9) is defined in prop. 4.1; 
the integral in q gives T'^' = -T"\ i 1, 2, 3; then T"' = 0, i = 1, 2, 3 
and this implies Toj = 0, i = 1, 2, 3. 

(ii) Set in H2.15|l a — i, [3 = j, i j; now compute the integral using the 
change of variables p — Af_i(0)g if {i = l.j ~ 2) or (i = l,j = 3) and 
p = Mi(f )g if {i = 2,j = 3), to obtain: T^' = -T'-?' if i ^ j; then 
yi2 ^ 2.13 ^ 7-23 ^ and = T13 = T23 = 0. 

(iii) Set in (|2.15|l a ~ (3 — 2 and compute the integral using the change of 
variables p — Mi{^)q. The integral in q gives T^^ — T^^ and hence 
^22 — 2I33. 

This completes the proof of proposition 4.1 ■ 

In all what follows, we assume that the collision kernel A satisfies assumption 
H4.6|l . and that f{t) is invariant by G . When we will study the coupled Einstein- 
Boltzmann system, we know by prop. 4.1 that it will be enough that /o be 
invariant by G. Also notice that the kernel A defined in remark 3.1 is an 
example of kernel satisfying assumption (|4.6|l . 

4.3 The Constraint Equations 

The Einstein equations (|2.14(l in which G = 1, give, using S"^: 

S° - 8ttT° = (4.10) 

It is proved in T, p. 39, that, in the general case, the four quantities: 

= — SttT" do not involve the second derivative with respect to t, of the 
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metric tensor ga/3 and, using the identities Va{S°'^ - SttT"^) = 0, that the 
quantities satisfy a linear homogeneous first order differential system. Con- 
sequently: 

1°) For t = 0, the quantities — 8ttT^ express uniquely in terms of the initial 
data ao, do, 60, ^0 and /q. 

2°) If (5*^ - 87rr^)(0) = 0, then (5° - 8TrT°)(t) = in the whole existence 
domain of the solution (a, b) of the Einstein equations l|4. 1(1 - 1)4. 2() - H4.3|l . In 
other words, ((4.10|l is satisfied if the initial data satisfied the constraint 

5°(0) = 87rT°(0) 

For this reason, the equations (|4.1U|) are called constraint equations. Notice 
that, since 5," = and T!^ = by prop. 4.1, equations H4.10|l for a = i, which are 
called momentum constraints, are automatically satisfied. It then remains 
equation (|4.1U|) for a = which is called the Hamiltonian constraint. Now 
5g - 87rTo° = is equivalent to: 

= 87:T°° (4.11) 

which is exactly H4.1|l . since 5oo = S"^"; Too = So, H4.1(l is the Hamiltonian 
constraint which is satisfied if it is the case for t — 0, i.e, if the initial data 
oq, 60, do, 60, /o satisfy, using H2.15|l for a = f3 — 0, 14. 111 . (|4.11|) . the initial 
constraint ^ 

=8vr/ p"{0)fo{p)aobldp (4.12) 

with, using ^ and /(O) = y/l + al{p^)^ + &g[(p2)2 + (p3)2]^ ap- 

pears that, the choice of the four independent initial data ao,bo,bo, fo uniquely 
determines do- We will take: 

ao > 0, 60 > 0, bo > 0, /o e Ll{R^), /o > a.e; do > 0. (4.13) 

where do > is obtained by taking bo sufficiently large. 

Notice that, in the 3+1 formulation of the Einstein equations, the Hamiltonian 
constraint 1(4. l|l writes: 

(trKf - K,jK'^ = 167rToo (4.14) 

where K = {Kij ) is the second fundamental form induced by the metric g on the 
hypersurfaces St — {t} x R^. Kij is given in the present case by Kij = —^dtgij, 
which gives, using ((2.11) 

Kn - -ad; K22 = K33 = -bb; K,, =0 if i^j (4.15) 

and trK — g''^ Kij is the trace of K; trK which represents the mean curvature 
of the space-time is then given by: 



trK = - 




(4.16) 
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4.4 The Evolution Equations 

The Einstein evolution equations are the second order non-Hnear o.d.e (|4.2|) - (|4.3|l 
in a and b. we study the initial values problem for this system on [0,T], T > 
with initial data defined by H4.4|l and satisfying the initial constraint H4.12|) . As 
we explained above, we are funded to consider the Hamiltonian constraint (|4.1|l 
or (|4.14() . as auxilliary equation. 

In order to apply standard results, we are going to show that the Einstein 
evolution equations are equivalent to a non-linear first order system. We con- 
sider the system (I4.1ll - H4.2|l - H4.r>|l and we set in the sources terms, with 
defined by (j^T^ : 

4.17) 

Pi = 4.18 

P2 = ^ (4.19) 

R = ^ ' ' ^ (4.20) 
P 

R+ = ^—^ (4.21) 
P 

Next we consider, following Rendall,A.D and Uggla,C. in ^5], the change of 
variables: 



3 
1 



52 + 2a2 
1 b 



""- = 111-' 
H is called the Hubble variable. Finally we set: 

P 



n = 



q^2Y.\ + ^{\ + R) 



4.22) 
4.23) 
4.24) 
4.25) 

4.26) 
4.27) 



17 is called the normalized energy density and q, the deceleration parameter. 
We have the following immediate consequences of the above definitions 
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Lemma 4.1 

< z < 1; < s < 1; 



2z 



0<Pi+2P2<p; 0<i?<l; 0<g<2 



r2 > 0; f7 = 1 - Sj 



(4.28) 
(4.29) 

(4.30) 



Proof 

• i) H4.23|) - (|4.24|l imply 0<z<l,0<s<l and solving this system in 
a^, 6^, yields the solutions given in (|4.28|) . Next, (/(t) > 0,a.e) implies 
p = 87rToo > 0; H4.26|l then shows that > 0. For the last result in 
consider (|4.14() that writes, using H4.17|l and (|4.26() : 

{trKf - KijK'^ = enH^, in which, by 



(4.31) 



Now use (EUSl), (O^ and ifi:^ which give f = -{trK)-2\, trK = -3H, 



and 



(E+ + 1)H to express f, | in l|OT|l and obtain = 1 - S 



• ii) To obtain the first resuh in use definition lfTT5|l 
of T"^, and (gSJ which gives T22 = T33. Next, use (|Onjl to obtain 
< i? < 1 and to obtain < g < 2 since > ^ < 1. 

• iii) Use expression H2.3|l of p^, expressions of a^, in (|4.28|l . the expres- 
sions of Pi, P2 given by H4.18(l . (j4.19|l to obtain (4.30). This completes 
the proof of Lemma 4.1 ■ 

Now the system (I4.1|l - (|4.2II - H4.3ll writes, using the notations H4.17|l - H4.18|l - 



ab 



(4.32) 



Pi 



(4.33) 



a 
a 



ah 
ab 



P2 



(4.34) 
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It shows useful to explicit the second derivatives. H4.34|l gives, using (|4.33|) to 



express |: 



a 
a 



-P2 



1 b 
2 2 U 



db 
ab 




2 db 1 db 

3 ab 3 ab 



(4.35) 



and we write | given by (4.33) on the form: 




2 / . \ 2 
Pi I I b\ lib 



3 \ b 



6 \ b 



(4.36) 



Now, if we use (|4.32|l to express the last terms in (I4.35f) and in (|4.36|) , we obtain 
the Einstein evolution equations in the form: 



db 
ab 



6 2^ ' 



(4.37) 




P^ 
2 



(4.38) 



We will also use the following relation, deduced from 14.16|l . (|4.22|l . H4.25|l : 



- = (l-2I]+)iJ 
a 



(4.39) 



Remark 4.1 Assume H > 0. Then: 

If bito) = at a given time t = to, then the problem of the global existence 
becomes trivial. 

The reason is the following: 
H4.36|l shows that 6 < 0, then 6 is a decreasing function. By (|4.28l) 

{n ^ I -Y,\>Q) ^ {-I <Y.+ < 1); hence, E+ + 1 > 0, and by b>0. 
But {b > 0,b decreasing, 6(to) = 0) ^ ( b{t) = for t > )• Hence b is con- 
stant and b{t) — 6(^0) for t > tg. 

Then, by l|051l . S+(t) = -1 for t > to, and by K7^ . n{t) = for t > to; 
then implies p{t) = 0, for t > to, and since < Pi + 2P2 < P, Pi{t) = P2{t) = 
for t > to- The Einstein equations l|4.32|l . (|4.33|l are trivially satisfied for t > to 
(since b = 0), and by d{t) = for t > to- Then: a{t) = Cit + C2; Ci, C2 

constants. 

Notice that {p{t) = for t > to) ^ (/(t) = for t > to). 

Conclusion: For t > to, the space-time becomes empty and the problem of global 
existence becomes trivial. 

In what follows, we will show that H > ii H{0) > 0, and we look for b 
such that 6 > 0; consequently, by (|4.25|l . E+ + 1 > 0. 
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Proposition 4.2 The Einstein evolution equations | |^.5'61 )- f^.5'7| ) in a, h, are 

equivalent to the following non-linear first order system in [H, s, z, 



i§ = -il + q)H^ 



(IV. 1) 



s)^+H 



dt 



6s(l 

:2z(l-z)(l + E+-3sI]+)iJ 



(4.40) 
(4.41) 
(4.42) 
(4.43) 



Proof: 



1. Suppose we have (|07jl - lll?^ . 



• i) and KT^ gi\ 



dH 
~dt 



1 i h b 



- 2 



which gives, using l|07jl - l|0^ . llOT|l and ^J^: 

dH I ,P Pi 1 r, ,T 

(4.40) then foUows from g2ni,(|i!23, and (|07|l . 

ii) Definition (4.24) of s gives 1 — s = ^^^^^^ and we can write: 



ds „ , ,6 
— = 2s(l - s - 



we then deduce from (4.25) and (4.39) that: 

ds , , f (l-2T.+)H 

-^2.(l-,s)(E, + l)H(^l-^^^-^ 

iii) Definition (4.23) of z gives: 



6s(l - s)T.+H 



dz r,2(^l ^ \ 

Tt^ \a^^ hh^j- 

(4.42) then foUows from (4.39), (4.25) and (4.28). 
iv) We can write, using (4.25) and (4.40): 



^-^d-tillbrb^'^'^^H 



(a) 
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Now the evolution equation (4.38) gives, using (4.25) and (4.39) 

h (bV b (id 4b\ p Pi _b p Pi 

b~\b) -b\i~a~Yb)-Q-^ -6^-1-2^+)^- 6 -T 

so that, (a) gives, using once more (4.25): 

— = -(,-2E+)- — - — = (S+ + l)(g-2E+)i7- — - — 
= (g - 2)I]+iI - lY.lE + qH-J^-^ 

Now express q in the third term by (4.27), use (4.26) that gives 
± = Mh;^ and definitions (4.20) and (4.21) of R and R+ to obtain 
(4.43). 

2. Suppose we have (4.40)-(4.41)-(4.42)-(4.43) . 

• i) Define the quantities a^, b^ as in (4.28) and the other quantities the 
same way as we did. Differentiating both sides of the corresponding 
formula (4.39) yields: 



= H{1 - 2S+) - 27JS+ 
which gives, using (4.40) and (4.41): 

- = -(1 + q){l - 2S+)if2 _ 2H^[{q - 2)S+ + 



Now, use expression (4.27) of q to obtain 



{l+R) + 2nH^R+ 
(b) 

For the second term of (b), use (4.26) which gives fiJf^ = ^ and the 
definitions (4.20) and (4.21) of R and R+ to obtain: 

-^(1 + i?) - 2nH^R+ = -| + i(P, _ 2P,) 

Next, for the first term of (b), use (4.25) to express S+, (4.16) and 
(4.22) that give if = + 2^) to obtain: 

2 F/'X^ 



ab 
ab 



and (4.37) follows. 
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ii) Differentiating both sides of | = (S+ + 1)H yields: 



which gives, using (4.40) and (4.43) and simpUfying: 



J = (-3S+ - 1)H^ - qH^ + QR+H^ + I ^ 



Now this writes using definition (4.27) of q: 

/ • \ 2" 



(-21]^ - 3S+ - 1)H'^ 



{l + R) + nH^R+ (c) 



For the second term of (c), use (4.26) which gives flH"^ = | and 
definitions (4.20) and (4.21) of R and R+ to obtain: 

flH 



Next, for the first term of (c), use (4.25) to express S+ and H = ^ + 2 
to obtain: 



(-21]^ - 3S+ - l)i?2 



db 
ah 



and (4.38) follows. This completes the proof of Proposition 4.2 ■ 

Now observe, using expressions (4.20), (4.21) and (4.27) of i?, i?+ and q, 
that these quantities express in terms of ^ and wo have nothing to bound such 
quantities. This is why, we give a new formulation of equations (4.40) and 
(4.43), that raise this problem, as follows: 



Proposition 4.3 Equations (4-40) and (4-43) can write respectively: 



Pi + 2P2 



(4.44) 
(4.45) 



Proof: 



• i) Write (4.40) using definition (4.27) of q in which R is given by (4.20); 
next, use (4.28) and (4.26) which give = 1 - and ^ = | to obtain 
(4.44). 
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• ii) Write (4.43) using definition (4.27) of q and apply definition (4.20) 
and (4.21) of R and next use once more (4.28) and (4.26) wfiicfi give 
]2 and ^ = 4. 

+ p 3H 



17 = 1 - S2 and ^ = to obtain (4.45) 



By Drop l4.3l system (IV. 1) is then equivalent to: 

^ + S^)ij2 _ (4 44) 

(IY2) I i=6.(l-.)S+i? (4.41) 

^ If = 2z(l-z)(l + E+-3sI]+)i/ (4.42) 
q± = -1(1 - S^)ifE+ + ii(E+ - 2) + ii(E+ + 1) . (4.45) 

It shows useful in what follows to introduce the following notation. 

Definition 4.1 LetO < s <1. Set: 

a(s) =/n/(s,l-s) (4.46) 

For the global existence theorem we will need the following a priori estima- 
tions. 

Proposition 4.4 Let 6 > be given and let to G M_|_. Then any solution 
(iJ, s,z,I]+) of the system (IV.2) on [ta,to + S], with H{to) > 0, satisfy the 
following a priori estimations, Vt G [0, S] : 



Hjto) 

i+H(to)/;:+*(i+g)dT 



^(^0 +t)= . . lit,. . . . (4-47) 



1 < ^.^e3H(to)t (4,48) 



Hito + t) - H{to) 

^ < , ] e'^(*»)* (4.49) 



a{s{to+t)) a{s{to)) 

^ < ] e'OHit'>)t (4 50) 



a{z{to + t)) Q;(z(io)) 

-l<S+(to+<)<^ (4.51) 

Proof: Let t £ [0,6]. By prop 14. 31 for H we consider (4.40) or (4.44), and 
for E+ we consider (4.43) or (4.45). 

1. (4.40) writes: —jj^ = {l + q)dt; integrating both sides on [t^, <o + A yields 
(4.47). Notice that, using (4.40) or (4.44), H \s a. decreasing function. 

2. By (4.29) we have < 9 < 2; so, by (4.47), (i/(to) > 0) ^ {H{to + t) > 0) 
and < H{to + t) < if (to)- Now (4.40) also writes: --^ = (1 + q)Hdt; 
integrating both sides on [to, to+t] and using < q < 2, H{to+t) < H{to) 
yields (4.48). 
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3. (4.41) gives: 



|6(1 - s)j:+H\dt; 



d(l - s) 



1 - s 



= \6sY.+H\dt 



Now integrate on [to, to +t], using |S+| <1, 0<s<l, 0<iJ< H{to) to 
obtain: 

r -6H{to)t <- s{to+t) ^ 6H(to)t 
I - s{to) - 

] -6H{to)t < l-s(fo+t) ^ g6fl"(to)* 
L — l-s(tn) — 

We then deduce, using definition (4.46) of a{s) that: 
a(s(to))e-s^(*")*<a(5(io + t)) 

and (4.49) follows. 

4. To obtain (4.50), proceed as above, using this time (4.42) which gives: 



= |2(l-z)(l+S+-3sS+)i7|dt; 



d(l - z) 



1 - z 



= \2z{l+T.+-3sY.+)H\dt 



Notice that the inequalities (4.48), (4.49) and (4.50) don't involve /, be- 
cause we used (4.40) in which < g < 2 and (4.41), (4.42), which contain 
no source terms of the Einstein equations. 



5. To prove (4.51), we show that, for every point ti G [to, to + 6], 



(S+(ii) = 2) 



i^+ih) < 0) 



(a) 



Since (a) will shows that i is an upper bound for all the values E_|_(t), 
t G [to, to + S]. We use equation (4.45) in S]_|_. Evaluating both sides of 
(4.45) Sit t ~ ti gives, using S+(ti) = |: 



(b) 



H' 16 4 2 

Now, the relation = 1 — gives fl{ti) — I — j = j. We then deduce 



from (4.26) that, lor t = ti, H'^ ^, and (b) gives: 



E+(ti 



1 

H 



-p-Pi+ 2P2 



(c) 



Now, the inequahties < Pi + 2P2 < P in (4.29) give : 2P2 < p - Pi so 
that (c) implies: 

S+(ti) = ;|(-y)<0 (c) 
since Pi > 0; and (a) follows. This completes the proof of proposition ^31 
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We deduce: 

Proposition 4.5 Let f £ C{Q,T; X^) , T > 0, be given. Suppose the initial 
value problem for the system (IV. 2), with initial data {H^, sq, zq,Ti+q) , with 
Hq > 0, < So < 1, < zq < 1, — 1 < E+o ^ ^, at t = 0, has a solu- 
tion (iJ, s, z, on [0,to] with <to < T. Then, any solution s,z,I]+) 
of the initial value problem (IV. 2) on [toi^o + S\, 5 > 0, with initial data 
{II,s,z,Yi+){tt)) — [H ,s,z,t]+){tQ), at t — to, satisfy, for t G [0,(5], the in- 
equalities: 

^ < 7oe^^°'^+') (4.52) 
< 7oe^^«(^+') (4.53) 



where 



H{to^ 


rt) ■ 


1 




a{s(to - 




1 




a{z(tQ 4 


-t)i 






70 = 1 





1 1 

So Zq 



(4.55) 



Proof: The hypothesis on / imphes that system (IV. 2) is defined on [0, T]. 
Apply prop 14. 4( Notice that by (4.40), iJ is a decreasing function; so (4.48), 
(4.49), (4.50) give in the present case, since -ff(to) < ^^o- 

{ 1 < _J_„3Hot. 1 < 1 ^6H„t 

H{to+t) - H(to)^ ' a(s(to+t)) - aisito))'^ jq^ ^^^^ 

a(z(to+t)) - aCzU))^^*^^"* 

Now apply prop l4.4l to the solution (7J, s, z,I]+) of (IV. 2) on [0,to], with 
(#,s,z,S+)(0) = (i?o,so,zo,i;+o); then (4.48), (4.49), (4.50) in which we take 
io = and t ^ to give: 

1 < J_g3Hoto. ^ < __}_^6Hoto. ^ <- _J_^WH„to ^|^) 



H(to) Ho a{s{to)) a{so) a{z{to)) a{zo) 

(4.52), (4.53), (4.54) then follow from (a), {h),to <T,t<5 and the definition 
(4.55) of 7o. This completes the proof of prop l4.5l B 

Corollary 4.1 Any solutions a, b of the Einstein equations (4.1)-(4-.2)-(4--3) 
on [to, to + 8\, to & (5 > 0, are increasing functions. 

Proof: (4.51) shows that 1 - 2E+ > on [to, to + S\; then (4.39) implies 
a > 0, since a > and H > Q; and a is an increasing function. Next, (4.25) 
gives: | = (11+ + 1)H and since 11+ + 1 > 0, H > 0, > this implies b > 
and b is an increasing function ■ 
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F{t, H, s, z, E+) 



Remark 4.2 By corollary \4-.l\ if we take in (4-4) ihe initial data oq, bo, such 
that 

3 3 
ao > 2? ^0-2 

then a and b satisfy the hypothesis of the existence theorem \3 . 1\ of solutions for 
the Boltzmann equation. 

We now prove the existence theorem for the Einstein equations, which are 
equivalent, using proD l4.3l to the system (IV. 2). Following Prop l4.4l we will 
take the initial datum Hq for H such that Ho > 0, and denote the initial data 
for the system (IV. 2) at i = 0, {Hq, sq, zq, S+q), where: 

Ho>0; < So < 1; < zo < 1; -1< S+o < ^ (4.56) 

We will apply the standard theory on the first order differential system. For this 
purpose, we will have to study the function defined by the r.h.s of the system 
(IV. 2), i.e, the function F defined by: 

6s(I - s)E+i7 
2z(l - z)(l + S+ - 3sS+)iJ 
-f (f - + ^(S+ - 2) + i|(S+ + 1) 

(4.57) 

(4.40) shows that iJ is a decreasing function, so that H{t) < Ho, t > 0, whenever 
H exists. We are then funded to suppose, using (4.28), (4.51), that F is defined 
for: 

(H, s, z, E+) e B =]0, Ho]x]0, l[x]0, l[x] - 1, i] (4.58) 

We will have to prove that F is continuous in t and locally Lipschitzian in 
X = {H, s, z, S+), with respect to the norm of we take to be 
\\X\\ = \H\ + \s\ + \z\ + \T,+ \. F depends on t through / in Pi and Pa (see their 
expressions in (4.30)); but, by hypothesis, / is continuous in t and so is F. 

A glance to (4.57) shows that the real problem will be to prove that Pi and 
P2 defined in (4.30) are locally Lipschitzian . 
We then begin by proving: 

Lemma 4.2 Suppose f e C {[to,to + 5]; Ll{R^)) , to>0, S>0. 
Let s, si,S2, z, zi,Z2 e]0, 1[, i?2 e]0. Ho], t e [to, to + S]. then: 

msi,zi,f ) - P.{s2,Z2,f)\ < ,6(,^),5(,,),8(,^),5(,^) ■ * = 1' ^ (4.59) 

a'^[s)a^(z) 



Pi{Sl,Zi,f) Pi{s2,Z2,f) 



Hi H2 
where C > is a constant. 



< g||/WII(|-'^l-^2| + kl-^2| + |gl-g2|) . ^ ^ ^ 2 

- H2{l + Hi)a^{si)a'^{s2)a^zi)a^Z2) ' 

(4.61) 
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Proof: Consider the expressions of p'^,Pi,P2 in (4.30). Let us begin by 
bounding the differences involved in Pi{sx,z\,f) — Pi{s2, Z2, f),i = 1,2. We 
have: 



1 



S2(l — Z2) Sl(l-Zl) 



1\2 



+ 



2Z2 



2zi 



(1-S2)(1-Z2) (1-Si)(l-Zi) 



p0(si,2:i)p0(s2,22) (pO(si,Zi) +P°(S2,Z2)) 



(a) 



Butp0(si,2i)p0(s2,^2)(p°(si,^i)+P°(s2,^2)) > (p°(si, ^i))V(s2, ^2) and (4.30) 
shows that (p^{s\,zi))'^p^{s2, Z2) is bounded from below by each one of the fol- 
lowing four quantities, that will serve to bound Pi and P2 respectively: 



si(l-zi) 1^82(1-22) 
2zi 



113 



(b) 



2Z2 



Sl(l-Zl) \(\-S2)(l-Z2) 

2zi ( 2Z2 



(T3ir)(T3i7) ( (i-JhVz2) ) ' b^l(b^)^ + (f^)^) 

quantities (b) will serve for Pi and (c) for P2. Next, we have: 



^ C'(|si-^2| + |zi-Z2|) 
— S1S2(1 — Zl)(l-22) 



S2(l — Z2) Sl(l — Zl) 
2Z2 



2zi 



(1-S2)(1-Z2) (1-Sl)(l-Zl) 



< 



C(|'il-S2| + |zi-Z2|) 
(1-Sl)(l-S2)(l-Zl)(l-Z2) 



(a) then gives, using (b)-(d) and (c)-(d) respectively: 



1 



1 



P°{si,zi) (52,2:2) 

1 1 



P°(si,Zl) P"(S2,Z2) 



< 



< 



C{\SI-S2\ + \ZI-Z2\) _l 
SlS2(l - Sl)(l - S2)2l2| 1^ 

C(|S1 - S2I + l-Jl - -22!) 1 



(c) 



(d) 



(f) 



SlS2(l-Sl)(l-S2)ziz| Ip'I 

We now bound the difFcrcnccs in the coefficients of the integrals in Pi , P2 which 
are given by (4.30). Concerning Pi, we have: 



1 



sf(l-,si) 
/ 1 



Zl 



1-21 



1 



S|(l-S2) 



Z2 



Z2 



,sf(l-Sl) S|(l-S2). 



1-21 



+ 



S|(l-S2) 



Zl 



l-zi 



Z2 



1-22 
(g) 
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in which: 



Sf(l-Sl) S|(l-S2) 



S|(l-S2)-Sf(l-Sl) 



(S2Sl)2(l - - S2) 



w4 - vs?)(i - S2) + sUsi - S2) 



(S2Si)2(l - Si)(l - S2) 



< 



C|si - S2I 



S?s|(l-Sl)(l-S2) 



and: 



Z2 



C 



(1^) - (1^) 



1-Z2 

5 



^1 



< 



1-Zi 
C\zi - Z2 



Z2 



1-Z2 



< 



zUi-Zl)^{l-Z2f 



SO that (g) gives: 



sf(l-si) VI -^1 



Zl 



S|(l-S2) 



Z2 



1-Z2 



< 



C{\Si - S2I + 1-21 - -221 



Sl{l - (1 - S2)Z^ (1 - Zi)^l - Z2f 



An analogous calculation gives, for the coefficient of the integral in P2 



Z\ 



Z2 



< 



C{\Si-S2\ + \Zi- Z2\) 

1 5 



Sl(l - Sl)2s|(l - S2)2^?(l - 2l)5(l - Z2)^ 

We can now establish (4.59). We have, using the expression of Pi in (4.30) 



(i) 



Pl{si,Zl,f) - Pl{S2,Z2,f) = 167r 



5^(1 -si) 



Zl 



l-zi 



1 



^2 



Notice that: 

p°{si,zi) > 



1\2 



S|(l-S2) 
1 1 



Z2 



Zl 



si(l - 21) 



P°{si,Zl) P'^{S2,Z2) 



\P'\ 



1-Z2, 

dp 



(j) 



(4.59) for Pi then follows from (h), (e), (j), \p^\ < ^/^TW, the definition 
(4.46) of a(s), < s < 1, and < a{s) < 1. 
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Now, concerning P2, we have, using its expression in (4.30): 

1 / Zl \^ 1 f Z2 



P2(si,Zl,/)-P2(s2,Z2,/) =32^ 



Notice that: 

''°"--'>( (i-„.T('i-..) )'i'''i « 



(4.59) for P2 then follows from (i), (k), \p^\ < y/T+W: (0 and (4.46). 

Next, (4.60) follows from the expressions of Pi, P2 in (4.30), using (j) for 
the integral in Pi and (k) for the integral in P2. 

Finally to obtain (4.61), we write, for i=l, 2: 

FT FT = TT " TT ■PH«i'^i'/) + 7r(^>*(*i'^i'-^)^^'(*2,22,/)) 

til tl2 V^l ^2 / ^2 

and (4.61) is a direct consequence of -g^ — = ^h^h^ ' (^•5'^)' (4-60) and 
< a(s) < 1 for < s < 1. This completes the proof of lemma I1.2B 

Now, concerning the other terms in (4.57), a usual calculation gives: 

f 1(1 + S^i)i/i2 - (1 + Y?^2)ill\ < CHoiHo + l)i\Hi -H2\ + - I]+2|) 

|si(l - Sl)S+iP-i - S2(l - S2)S+2i/2| < C{\Hi - H2\ + \si - S2I + |S + i - E+2I) 
|zi(l - Zi)(l + E+1 - 3siS + i)i7i - Z2(l - Z2)(l + S+2 - 3s2S+2)-ff2| < 

C{\Hi - H2\ + |si - S2I + \zi - Z2\ + |S+i - S+2I) 

1(1 - S^i)i/iS+i - (1 - S^2)^2S+2| < C{1 + H^){\Hi - H2\ + |S+i - S+2I) 

^Vsi, 32,21,22 e]0, i[;i/i,i/2 e]0, Ho];S+i,l]+2 e] - 1,5]- 

(4.62) 

We then deduce from (4.59), (4.61) and (4.62), using the expression (4.57) of F 
the inequality: 

||P(t,i7l,Si,2i,I] + i)-P(i,ff2,S2,Z2,S+2)|l < M(|i7i-H2| + |si-.52| + |2l-22| + |S + l-S 

(4.63) 

where: 



M = M(i7i,i/2,Sl,S2,2l,22) = C 



H2{1 - Hi) J a^{si)a^{s2)a^(zi)a^Z2) 

(4.64) 

We can now state: 

Proposition 4.6 Let f e C{[0,T]; Xr), T > 0, be given. Then, the initial 
value problem for the system (IV. 2), with initial data {Hq, sq, zq, S+o) at t — 
satisfying (4.56) has a unique solution {H, s, z,Y,^) on [0, P]. 
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Then: 



and: 



Proof: Take, following the definition (4.58) of the domain of F: 

(i?°,s",2°) e]0,i7o]x]0,l[x]0,l[ (4.65) 

and: 

s° + 1 z° z" + 1 

H. e] — ,Ho]; e]-, ^[; e]-, ^[; t ^ 1,2 (4.66) 

s" s" + l zO zO + 1 

H.>—;-<s.< y < < ^ = 1,2 

Consequently, the number M defined by (4.64) is bounded on the neighborhood 
] ^ , ffo] X ] 4 , ^ [ X ] ^ , ^ [ of (7? , 5° , zO ) in by a number M° , depending 
only on H°, s°, z°, r, since \\f{t)\\ < r, Vi € [0,T]; the inequality (4.63) then 
shows that F is locally Lipschitzian with respect to the norm of M.*. By the 
standard existence theorem for the first order differential system, the initial 
values problem for system (IV. 2) has a unique local solution {H, s, z, S+). Now 
proposition ^31 in which we take io = 0, t G [0, T], shows, using definition (4.46) 
of a{s), < s < 1, that the functions jj, i, jz^, -j, jz^ are uniformly bounded 
and (4.60) implies that the functions Pi, P2 are uniformly bounded; so F defined 
by (4.57) is uniformly bounded and, by the standard theory on the first order 
differential system, the solution {H, s, z,Y^^) is global on [0,r]. This ends the 
proof of DroD l4.6l B 

It shows useful to deduce the following result: 

Proposition 4.7 Let f e C'{[0,T]; X-r), T > 0, be given. Let (ff,s,z,S+) be 
the solution of the initial value problem for the system (LV.2) with the initial 
data {Hq, So, zq, S+q) att = 0, satisfying (4-56). Let to € [0, T[. Then the initial 
value problem for the system (IV. 2), with initial data {H, s, z, E+)(io) att — to, 
has a unique solution [H, s, z, E_|_) on [tg, io + ^\ where S > is independent of 
to- 

Proof: The proof of the existence of a local solution {H, s, z,Y^^) on an 
interval [to, to + S], for the initial value problem for the system (IV. 2) with 
initial data {H, s, z, T,+ ){to) is analogous to that of prop 14. 61 Suppose we look 
for 5 such that < (5 < 1; then applying the inequalities of prop 14.51 leads to an 
inequality (4.63) with a constant M independent of the initial data ai t = to, 
and of to', from there, the existence of such a number 6 > U 

We end this section by the following result: 

Theorem 4.1 Let f G C{[0,T]; Xr), T > 0, be given. Then the initial value 
problem for the Einstein equations (4.1)-(4-2)-(4-3) with initial data 
{ao , bo , CLo , bo , fo) att~0, satisfying (4.13), and the initial constraint (4.12), 
has a unique solution a, b on [0, T] . 
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Proof: Apply proposition 14. 61 choosing by virtue of the change of variables 
(4.22), (4.23), (4.24), (4.25) and using (4.16), the initial data iJo, sq, zq, S+o, 
as follows: 

(4.68) 



ap I 2 Op 
3 V 10 ba 



ag ^+26,/+!' 



V _ 1 6o 1 



Notice that ^ 

ao 



> 



SO 



^ ^ < |i?o; then -1 < E+o = - 1 < i 

that {Hq, So, zqj'S+o) defined by (4.68) satisfy the assumption (4.56). Proposi- 
tion ^21 then proves the existence of the unique solution {H, s, z, E+) on [0, T], 
of the initial value problem for the system (IV.2) with the above initial data 
(iJo) So, ^0, S+o) at i = 0. The equivalence of the Einstein evolution equations 
(4.37)-(4.38) with the system (IV.2), then shows that, (a, 6) defined in (4.28) 
in lemma 14.11 is the unique solution on [0, T] of the initial value problem for 
the Einstein evolution equations. Now, since the initial constraint (4.12) is sat- 
isfied, this implies that the Hamiltonian constraint (4.1) is satisfied on [0,r]. 
This ends the proof of theorem 14.11 ■ 



5 Local Existence theorem for the Coupled Einstein- 
Boltzmann System 

As we proved in paragraph 3 and paragraph 4, the Einstcin-Boltzmann Sys- 
tem in (a, b, f) is equivalent to the following first order differential system in 
(/, H, s, 2, E+) given by (3.15) and (IV.2): 



(V) 



dH _ _3/-| 
dt 2V-^ 

I = 6,s(l - s)j:+h 

§^2z(l-z)(l 



P1+2P2 

6 



T.+ - 3sS+)i/ 



dt 



(1-S2)i/E^ 



(5.1) 
(5.2) 
(5.3) 
(5.4) 

2) + ii(I]+ + l). (5.5) 



in which Q is the collision operator defined by (3.1)-(3.2)-(3.3) and Pi, P2 
are defined in terms of s,z,f by (4.30). To solve system (V), we apply the 
standard theory of the first order differential system for functions on R, with 
values in the Banach space E = Li,{R^) xMxRxRxK, whose norm will be 
taken to be: \\{f, H, s, z,1:+)\\e = \\f\\ + \H\ + \s\ + \z\ + |E+|. Since / is now 
also an unknown, the problem will be to prove that the function F, defined by 
the r.h.s of (V) i.e: 



Fif,H, 



-1(1 + 
6s(l- s)E+ 
2z(l~z)(l 



P1+2P2 

6 



H 

f S I - 3sE 



-i(l-I]2.)i/E+ + ii(E^ 



1) 



(5.6) 
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which does not depend exphcitly on t, is locahy Lipschitzian in (/, H, s, z, 
with respect to the above norm of E. FoUowing (4.58), F will be defined on: 



L^(K3)x]0,iJo]x]0,l[x]0,l[x] 



1, i]. A glance to (5.6) shows that, taking into 



account the study we did in paragraph 4 for the system (IV.2), the new problem 
to face here will be the study of all the differences in /, and the differences of 
^^^u^^ in s and 2, which turns out to be the only real problem we meet now, 
since expressions Pi, P2 in (4.30) show that these functions arc linear in /. We 
then begin by proving: 

Lemma 5.1 Let /i, /2 G L2{M.^), si, S2, zi, Z2 G]0, 1[ then: 



1 



1 



p"(si,Zl) P"{S2,Z2) 
1 



Q(/l,/l)(s2,22) 



< 



(Q(/l,/l)(s2,^2)-Q(/2,/2)(s2,^2)) 



C\\fin\si-S2\ + \ZI-Z2\) 

a^{si)a'^{s2)a'^{zi)a^{z2) 

(5.7) 

C{\\h\\+\\f2\\)\\fl-f2\\ 



P°{S2,Z2) 

where C > is a constant. 



< 



(52)0^(2:2) 



(5.8) 



Proof: For (5.7), we follow the proof of (4.59) in lemma using this time 
in that proof (/(si,zi))2 > and (/(si, zi))^ > ^^_^L_^[(p2)2+ 

{p'^)^] and, the relations (a) and (d) to obtain: 



P°isi,zi) 



P°is2,Z2) 



< 



C(|Sl-S2| + |2l-2;2|) 



P°{s2,Z2)zi{l~ Zi){l- Z2) \S1S2 (1-Si)(l-S2) 

(a) 

Now, we use (3.27) with / = /i, .g 0, (3.28) and (4.28) which gives the 
expression of a, b in terms of s and z to obtain: 



QiflJl)is2,Z2) 



P"(S2,Z2) 



< 



C\\fi 



S|(l-S2)(l-Z2)^ 



(b) 



(5.7) then follows from (a) and (b) and definition (4.46) of a{s). Now (5.8) is 
given directly by (3.27) with f = fi, g = f2, the expression of ab^ in terms of 
s, z, and once more (4.46). This completes the proof of lemma I^TTl M 
We will also need the following result which is not obvious: 

Lemma 5.2 Let f e L2{^^), si, S2, zi, Z2 g]0, then: 



1 



(Q(/,/)(si,Zl)-Q(/,/)(s2,Z2)) 



P°{si,zi 
where C > is a constant. 



< 



C\\fn\si-S2\ + \ZI-Z2\) 

a^{si)a^{s2)a'^{zi)a^{z2) 

(5.9) 



Proof: Here, we split Q into Q using the basic definition (3.1) of Q, 
i.e Q = — Q~ where Q"*", Q~ are defined by (3.2) and (3.3) in which ab"^ is 
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expressed in terms of s and z. It is here that we need the Lipschitz continuity 
assumptions (3.7), (3.8) on the coUision kernel A. We write: 

Q{f, f){si,Zl) - Q{f, f){S2, Z2) = [g+(/, /)(S1, Zl) - Q+(/, /)(.S2, Z2)] + 

[Q-{f,f){s2,Z2)-Q-{f,f){s^,z^)] (a) 
We can write, using the expression (3.2) of Q+i 



Q+(/,/)(Si,Zi)-Q+(/,/)(s2,Z2) = 



Q.{si, S2,zi,Z2,p, q,p', q')f{p')f{q')dqduj 



where f2 is a function of the indicated arguments, we write on the form: 

" = {S1,ZI)[A{SUZI)-A{S2, Z2)]+ (Sl, Zl) - {S2, Z2) 



(b) 



A{S2,Z2) 

(c) 



in which A{si,Zi) stands in fact for A{a{si,Zi),b{si,Zi),p,q,p',q') i=l, 2; we 
now boimd the first term in (c). The assumptions (3.7), (3.8) on the colhsion 
kernel A give, adding and subtracting A{a{s2, Z2), 6(si, ^i)): 

\A{si,zi) - A{s2, Z2)\ < A;o(|a(si, zi) - a{s2, ^2)! + \b{si, zi) - b{s2, Z2)\) (d) 

Now (4.28) gives, by usual calculation: 



\a{si,zi) - a{s2,Z2)\ 



\b{si,zi) - b{s2,Z2)\ = ^J\ 



2zi 



2Z2 



(1-Sl)(l-Zl) V (l-«2)(l-22) 

We then deduce from (d), (e) and (4.28) which gives: 

3 

2 



^ C(l«l-^2l + kl-Z2l) 
^ 1 

8iS2zl (l-2l)(l-Z2) 

C(|S1-^2| + |Z1-Z2|) 



< 



{ab'){s,z) 



S2 



1 - z 



(1-Sl)(l-S2)2f (l-Zl)(l-22) 

(e) 



(f) 



that: 



-g- ] {si,Zi)[A{si,Z-i) - A{S2, Z2)\ 



< 



C{\si - S2I + \zi - Z2I 



a^{si)a^{s2)a^{zi)a{z2)q°{si,zi) 



Now, we bound the second term in (c); we write: 



(g) 



ab' 



(si,zi) 



ab' 



(S2,Z2) 



(a6^)(si,0i)g°(s2,Z2) ^ {ab'^){s2,Z2)q'^{si,Zi) 
q°{si,Zi)q"{s2,Z2) 

{ab'){suz,)- '^yyl'l {ab'){s2,z2) 



q°{si,zi) 



9°(si,^i) 



{ab'^){si,zi)-{ab'^){s2,Z2)+ 1- 



9°(S2,22) 

q°{si,zi) 

q°{S2,Z2) 



{ab'^){s2, Z2) 
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Hence: 



q°{s2,Z2) - (?°(.Si,Zi) 



Now, by (4.30): 

9°(-S2,^2) - g°(si,Zl) _ ( 



82(1-22) si(l-zi) 



1\2 



g°(s2, ^2) 9°(S2, 22) [(7°(S1, Zl) + g°(s2, Z2)] 



+ ' gO(s2,z2)[gO(si,zi) + gO(s2,Z2)] 
and we proceed as in the proof of (4.59), using this time: 



[q'is2,Z2)Y> 



yr-, [q°iS2,Z2)Y> 



82(1 - Z2) 

to obtain, using expression (f) of afe^: 

?°(S2,22) - g°(si,^l) 



(a6 ) (52,2^2) 



9° (52,22) 



< 



(l-S2)(l-22) 

C(|S1-.S2| + |Z1-Z2|) 

a2(si)a3(s2)Q;2(0i)a3(2;2) 



(i) 



Now, concerning the first term in (h), we have, using (f) and by usual calculation: 



|(a62)(si,^i)-(a62)(s2,22)| =2 



< 



C(|si - S2I + |zi - Z2I 



a (s 1 ) a (.S2 ) a''' (zi )q:3 (z2 ) 
We then have for the second term in (c), using (h), (i), (j) and < A < Cq: 



(j) 



^(S2, Z2) 



a}? 



(Sl,2l)-(^^)(S2,Z2) 



< 



C(|si-S2| + kl-22|) 



(j°(si, zi)a2(s^)Q;3(|g2)a^(zi)a3(2;2) 

(k) 

We then deduce from (b), using (c), (g) and (k) that: 



{Q+{.fJ){si,Zl)-Q+{f,f){s2,Z2)) 



P"{si, Zl) 

C{\SI^ S2\ + \ZI^ Z2\) 

a^{si)a^{s2)a'^{zi)a^{z2) 



< 



VT+W fiP')\\f{Q')\ 



S2 p°q°{si,zi) 



dpdqduj (1) 
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We compute this integral the same way as in the proof of (3.25) in lemma ISTI 
using the change of variables {p, q) i-^ {p', q') defined by (3.12) to obtain: 

q|/|P(|si-S2| + |zi-Z2|) 



1 



(Q+(/,/)(si,Zi)-Q+(/,/)(,S2,Z2)) 



< 



a^{si)a^{s2)a^{zi)a^{z2) 
(m) 

We now proceed the same way for the second term in (a), i.e 
{Q~{L f){s2,Z2) — Q~{f, f){si, zi)); the only difference is that, in the integrals 
(b) and (1), f{p')f{q') is replaced by f{p)f{q)- So, no need this time to change 
the variables and a direct calculation using p'^q'^ > 1 leads to the same estimation 
(m), just substituting to (5+, and lemma \S7^ follows ■ 

Now we deduce from lemma ISTI and lemma l5^ the following: 

Lemma 5.3 Let fi, /2 G Ll(R^), si, S2, zi, 22 g]0, 1[; then: 



0(/l,/l)(Si,Zi) g(/2,/2)(s2,;^2) 



P°(si,2i) 



P°is2,Z2) 



< 



C^(ll/l|P + ll/l|| + ll/2|l) 



a^{si)a^{s2)a^{zi)a^{z2) 

(||/l-/2|| + |si-S2| + kl-^2| 



(5.10) 



where C > is a constant. 
Proof: Write: 

Q(/l,/l)(si,Zl) Q(/2,/2)(s2,Z2) 



QifiJi)isi,zi) - Q(/i,/i) (52,^2) 



P°{si,zi) 



P°{S2,Z2) 



P°{si, Zl) 



1 



1 



p"(si,Zi) p'-'[S2,Z2)J 

^ Q{hJl){s2,Z2)^Q{f2j2){s2,Z2) 
P^{S2,Z2) 

and apply: (5.9) with / = /i to the first term, (5.7) to the second term, and 
(5.8) to the third term and < a{s) < 1, to obtain (5.10) ■ 
Finally, concerning the differences in Pi and we prove: 

Lemma 5.4 Let fi, f2 G Ll{R^); si, S2, zi, Z2 e]0, 1[; Hi, H2 e]0,i?o]- Then: 

C(l + ||/2||)(||/l-/2|| + kl-S2| + kl-Z2|) 



mSuZiJi)~P,is2,Z2j2)\ < 



-R(si,Zl,/l) fj(s2,Z2,/2) 



Hi 



H2 



< 



a^{si)a^{s2)a^{zi)a^{z2) 

(5.11) 

C{1 + ll/lll) (ll/l - /2II + \SI - S2I + \ZI - Z2\ + \Hi - H2\ 



1=1, 2; where C > is a constant. 



H2(l + Hi)a^{si)a^{s2)a^{zi)a^{z2) 

(5.12) 



Proof: The expressions of P^, i=l, 2 in (4.30) show that Pi is linear in /. 
So we can write: 

Pi{si,Zi,Ji)-Pi{s2,Z2,j2) = Pi{si,Zi,Ji-j2) + {Pi{si,Zi,j2) - Pi (S2 , ^2 , /2 ) ) 

(a) 
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Pi{Sl,ZiJi) Pt(s2,^2,/2) ^ _^ / Pi{si,ZiJi) _ -Pi(s2,Z2,/l) \ ^ Pt{s2, Z2J1 - 

Hi H2 \ Hi H2 ) H2 

(b) 

so, to obtain (5.11), apply (4.60) with /(t) ~ fi — f2 to the first term of (a), and 
(4.59) with f{t) = /2 to the second term of (a). To obtain (5.12), apply (4.61) 
with f{t) — fi to obtain the first term of (b), and (4.60) with f{t) = /i — /2 to 
the second term of (b) ■ 

We can then write, using (4.62), (5.10), (5.11), (5.12) and the expression 
(5.6) of ^: 

||F(/i,i/i,Si,Zi,S + i) - F{f2,H2,S2,Z2,^+2)\\E < 

N (ll/i - /2II + \si -S2\ + \zi- Z2I + \Hi - H2\ + - E+2I) (5.13) 
where: 

N = Nifi,f2,Hi,H2, si, S2, zi, Z2) 

r \, I f, I 1 \ A1 + II/1II + II/2II + II/1IP) 

[ \ H2{\ + Hi)) \a^{si)a^{s2)a^{zi)a^Z2) 
We can now state: 

Proposition 5.1 There exists a number I > 0, such that, the initial value 
problem for the differential system (V), with initial data (/q, Hq, sq, zq, S-|_o) at 
t = {), satisfying (4-56) and /o G L2(R^), has a unique solution {f , H, s, z,Yi^) 
on [0,1]. 

Proof: Choose H° , 5°, z° as in (4.65) and /" e L\{W^). 
Take Hi, Si, zi, i=l, 2 satisfying (4.66) and 

/i, /2 e B{f, 1) = {/ e Li(M3), 11/ - ni < 1}. Then H„ s„ z„ i=l, 2 satisfy 
(4.67) and ||/,|| < ||/"|| + 1, i=l, 2. Consequently, the number N defined by 
(5.14) is bounded in the neighborhood B{f° ,l)x]i^, Ho]x]^, ^[x]^, ^[ 
of (/° , ° , s° , z° ) in (R3 ) X R3 , by a number N° , depending only on f° , H° , s° , z' 
The inequality (5.13) then shows that F is locally Lipschitzian with respect to 
the norm oi E = Ll{M.^) x R'* and proposition |0] follows from the standard ex- 
istence theorem for the first order differential system, for functions with values 
in a Banach space. Notice that / G C ([0, ^2(1^^)) ■ 
We end this paragraph by the following result: 

Theorem 5.1 Let fo £ Ll{M.^); gq, bo, oq, bo satisfying (4.13), the initial 
constraint (4.12), and oq > §, ^ §> be given. Let r > ||/o||. 
Then, there exists a number I > 0, such that, the initial value problem for 
the coupled Einstein- Boltzmann system (2.13)-(4.1)-(4-2)-(4-3) has a unique 
solution /, a, b on [0, 1]. The solution (/, a, b) has the following properties: 

a and b are increasing functions (5.15) 

f eC{[0,l];Xr) (5.16) 

Ill/Ill <ll/o|| (5.17) 
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Proof: Apply proposition l5.ll choosing /o £ Ll{R^), /o > a.e, and define 
Ho, So, zo, S+o by (4.68); (4.56) is then satisfied. The existence of a unique 
solution {f,a,b) of (2.13)-(4.1)-(4.2)-(4.3) on [0,/], Z > 0, is a direct consequence 
of proposition l5.1l the equivalence of the above system and the system (V), and 
the hypothesis on the initial data at t = 0. 

Now, concerning the properties of the solution (/, a, b): 

• i) (5.15) is given by coroUarv 14.11 with to = 0, S = I. 

• ii) The hypothesis ao > |, 6o > | and (5.15) show that a{t) > |, 
b{t) > |, Vt e [0, a and b then satisfy the hypotheses of the exis- 
tence theorem 13.11 of the Boltzmann equation in /, in which we take 
to = and T — I > 0. (5.16) is then a consequence of the uniqueness and 
ofC([0,Z];X,)cC([0,/];Li(R3)) 

• iii) (5.17) is given by (3.19) with to = ■ 

6 The Global Existence Theorem for the Cou- 
pled Einstein-Boltzmann System 

6.1 The Method 

We assume, in all what follows, that the initial data (/o, ap, bo, ao, &o) at t = 0, 
for the Einstein-Boltzmann system (2.13)-(4.1)-(4.2)-(4.3) satisfy (4.12), (4.13), 
o-Q ^ ^, bo > |, and that, in the initial data (/o, Ho, sq, zo, S+o) at t = 0, for the 
system (V), paragraph 5, Ho, Sq, zo, S+o are defined by (4.68), and hence satisfy 
(4.56). Denote [0,T[, T > 0, the maximal existence domain of the solution, 
denoted here (/, _ff , s, z, E+) and given by proposition 15.11 of the initial value 
problem for (V), with the above initial data at t = 0. 

If T = +00, the problem of the global existence is solved. We are going 
to show that, if we suppose T < +oo, then the solution (/, H, s, z, can be 
extended beyond T , which contradicts the maximality of T. 

The strategy is the following: 
Suppose < T < -l-oo and let to E [0,T[. We will show that there exists a 
strictly positive number 6 > 0, independent of to, such that, the initial value 
problem for the system (V), on [to, to + S], with initial data 
if (to), H(to), s{to), z{to), S+(io)) at t = to, has a unique solution (/, H, s, z, I]_|_) 
on [to, to+i^]- Then, by taking to sufficiently close to T , for example, to such that, 
< T - to < f , hence T < to + |, we can extend the solution (/, H, s, z, S+) 
to [0,to -l- f [, which contains strictly [0, T[, and this contradicts the maximality 
of T. In order to simplify the notations, it will be enough to look for a number 
6 such that < (5 < 1. As in previous paragraphs, we denote r, a real number 
such that r > 1 1 /o 1 1 • 
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6.2 The functional framework 

Proposition 6.1 Let G [0, T[ and < S < 1. Then, any solution (/, H, s, z, E4_) 
for the initial value problem for the system (V) on [ioi^o + with the initial 
data at t — t^: 

if, H, s, z, I]+)(io) = (/>o), i?(to), .5(to), 5(io), S+(io)) (6.1) 
satisfy the inequalities: 

7777^—7 < Mo; , < ^^o; , ^ < ^o; t e [0, S] (6.2) 

H{to + tj a[s{to+tjj a(z[to + t)) 

where: 

Mo = Mo(ao, bo, ao, bo, T) = ( ^ + - + -) e^'^'oiT+i) (g.S) 

V Ho So zo J 

in which Ho, so, zo are defined in terms of ao,bo, ao,bo by (4-68). 

Proof: Apply proposition ^3] to the subsystem (IV. 2) of (V), and consider 
the solution {H, s, z, S+) = {H, s, z, S+) of the initial value problem (IV.2) on 
[0, T[, with initial data [Ho, So, zo, S+o), defined by (4.68) in terms ag, ho, ho, bo', 
since to e [0,T[ and {H,s,z,^+){to) = [H ,s,z,t+){to) = {H ,s,z,t+){to), 
(6.2) and (6.3) are given by (4.52), (4.53), (4.54), (4.55), (4.68), and < (5 < 1 
■ 

In all what follows, Mq is the absolute constant defined by (6.3). We deduce 
from (4.46), (6.2) and the expressions of a^, 6^ in (4.28) that: 

- < Mo; - < Mo; < Mh b^ < 2M^ (6.4) 

s z 

Wc then deduce from (6.4), the definition (4.23), (4.24) of z and s, in terms of 
a and b, using > Cq > | > 2, and the inequality for H in (6.2): 

1 11 1 11.... 

< z < TT-, < s = ^ — ^ < Tj; < H < Ho 

Mq- - 1 + 2Mo^ Mo - 1 + 2a25-2 - i + 2Mo^ Mo " " 

(6.5) 

On the basis of (6.2), (6.3), (6.5) we now introduce the following functions 
spaces, for to G [0,T[ and S > 0: 

K = e C{to, to + 5);^<s< I ; / < Mq; t e [0, S] 

[ Mo l + 2Mo a(s{to+t)) 



= \He C{to, to + 5);^< H{to + t)<Ho;te [0, 5] \ (6.6) 



1 

Mo 

Gl = |e+ e C{to, to + 5); -1 < S+(io + *) < i; te [0, 5] 
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One verifies easily that E^^^, i^j?, Gf^ are complete metric subspaces of the 
Banach space denoted C{to, to + S), of the continuous (and hence bounded) 
functions on the line segment [to, to + S], endowed with the norm: 

\H^= sup \uit)\; ueCito,to + 6) (6.7) 

tG[to,to+<5] 

6.3 The global existence theorem 
Proposition 6.2 Let to £ [0,T[. 

There exists a strictly positive number 5 > 0, depending only on the absolute 
constants ao,bo,do,bo,T and r, such that, the initial value problem for the sys- 
tem (V), with the initial data (/(to)i H{to), s(io), ^(^0), S+(to)) at t — to, has a 
unique solution (/, H, s, z, S+) € C([to, to + S]; Xr) x x x Ef^ x G^. 

Proof: By theorem 13.11 we know that, if we fix s, z G Ef^^ and if we define 
a = a{s, z),b — 6(s, z) by (4.28), then equation (5.1) in / has a unique solution 
/ e C{[tQ, to + S];Xr) such that f{to) = /(to), and, by (3.19) and (5.17) with 
/ = /, that: 

||/(t)|| < ||/>o)ll < ll/oll <r, te[to,to + S] (6.8) 

Next, by proposition 14.71 we know that, there exists a number S > 0, (we 
can suppose < S < 1), such that, if / is given in C{[to, to + S]-Xr), then 
the system (IV. 2) has a unique solution (i7, s,z,I]+) on [to, to + S] such that 
{H, s, z, E_|_)(to) = {H, s, z, S+)(io). Now Droposition l6 . II and inequalities (4.51) 
show that {H,s,z,T.+) e xF^^ x Ef^ x Ef^ x Gf^. This allows us, setting 
Xl = G{[to, to + S];Xr) X {e\ x EI^) and 

Y^l = G{[tQ, to + S];Xr) X {Fl x eI^ x Ef^ x GfJ, to define the appfication: 

g : Xl -> r/j (7, {s,z)) ^ (/, {H, s, z, E+)) (6.9) 

We are going to show that, we can find i5 > 0, such that g defined by (6.9) 
induces a contracting map of the complete metric space Xf^ (defined above) 
into itself that will hence have a unique fixed point {f,s,z). This will allow 
us to find H and S+ such that (/, (i?, s, z, E+)) be the unique solution in Y^^ 
(defined above) of the initial value problem for the system (V) with the initial 
data (6.1). 

So, if we evaluate the r.h.s of (5.1) for, s = s, z = z £ Ef , and if we set, in 
(5.3) and (5.5) P, = P, = P,(s,zJ), i=l, 2; where / € C{[to\ to + 5];Xr), then 
there exists a solution (/, (7J, s, z, E+)) of that system, taking the initial data 
(6.1) at t = to, or equivalently, a solution of the following integral system: 



(VI) 



f /(^o +t) = /(^to) 

H{to + t)^H{to)+ 



■Qif, f)is, z){T)dT 
(l + E^)i/2 



P1+2P2 

6 



(T)dr 



s{to + t)^ I{to) + /;°'*' 6s(l - s)Y.+H{T)dT 



z{to + t)^z{to)+ //;+*2z(l-z)(l + 



3sE+)iJ(T)dr 



^+{to + t) = j:+{to) + -|(i-E^)i/E 



i^(s+-2) 



||(E+ + 1) (r)dr 



40 



To (/i, {si,Zi)) G Xf^, i=l, 2, corresponds the solution 
(/j, {Hi, Si, Zi, S+i)) G i^Q, i=l, 2, of the above integral system (VI). Writing the 
equations for i=l, i=2 and subtracting yields, using respectively: (5.10) with 
s, = .s„ = z„ IlLMll < r; (5.11) and (5.12) with fi = h \\\fi\\\ < r; (4.62), 
definition (6.6) of Ef^, F^, (6.7), < t < S: 

lll/i - Mil < SM, [|||/i - Mil + ||5i - 52|U + 11^1 -^2|U] (6.10) 



m-H4^<6M2\\\f,-f,\\\ 

+ SM2 [||i?l - H2\\^ + - S2|L + 11^1 - ^2|L + - S+2IIJ 

(6.11) 

||S1 - S2IL < SMs [\\H, - F2IL + \\S1 - S2II00 + 11^1 - ^2||oo + IIS+1 - S+2IU] 

(6.12) 

11^1 - ^2|L < <5M4 [||i?l - H2\\^ + \\S, - S2IL + Ikl - ^2|L + l|5^+l - S+2IL] 

(6.13) 

||S + l-S+2||^<5M5|||7i-72||| 

+ 5Ms [\\H^ - if2|U + - S2IU + Iki - ^2|i^ + llS+1 - E+2II J . 

(6.14) 

where Mi, M2, M3, M4, M5 are absolute constants, depending only on the 
absolute constants ap, bo, ciq, bo, T and r. 

We keep (6.10) and add (6.11)-(6.12)-(6.13)-(6.14) to obtain: 

[\\H, - i/2lL + ll-'^i - 52IL + Iki - ^2|L + 11^+1 - S+2IIJ 

<4(5(M2 + M3 + M4 + Ms) II |7i - 72 1 II 

+ 4<5(M2 + M3 + M4 + Ms) [||ifi - + - S2II00 + INi - ^2|L + l|5^+i - S+2IIJ . 

(6.15) 

If we choose 5 such that: 

16(Mi+M2 + k + M4 + Ms) ) ^'-''^ 

which imphes that 5Mi < \ and 4^(M2 + M3 + M4 + M5) < \; then (6.10) and 
(6.15) give: 

/ lll/l - /2III < ^ - S2IL + -~2||oo] 

\ [||Fi - H^w^ + ||si - S2IU + 11^1 - ^2|L + lis+i - S+2II J < 5 II l/i - /2I II • 

and by addition: 

lll/l - /2III + [\\H, - F2IL + - S2IL + INl - -^2110, + 11^+1 - ^+2|l J 

<^[|||7l-72|||+l|5l-S2|L + ||zi-Z2||J 
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from which we deduce: 

lll/l - MII + IISI - S2|L + ||Z1 - Z2IL < I [|||7l -72III + 11^1-5211^ + \\Z1-Z2\U 

(6.17) 

(6.17) shows that, the map (/, (s, z)) (/, (s, z)) is a contracting map from 
the complete metric space Xf^^ into itself, for every number S satisfying (6.16), 
which shows that 6 depends only on the absolute constants Mi, i=l, 2, 3, 4, 5. 
This maps has a unique fixed point (/, (s, z)) € ^to- Now, to determine H and 
E_|_, since s is known, (5.3) determines the product in HT,^ in terms of s; then, 
substituting this product in (5.4) gives H in terms of s and z; once H is known, 
(5.3) gives E+, and we obtain the desired unique solution (/, {H, s, z, S+)) e Y^^. 
This ends the proof of proposition 16. 21 B 

We can then state: 

Theorem 6.1 The initial values problem for the spatially homogeneous Einstein- 
Boltzmann system on a locally rotationally symmetric Bianchi type I space-time, 
has a global solution (/, a, b) on [0, +00 [, for suitable arbitrarily large initial data 
att^O. 

Remark 6.1 1. Nowhere in the proof we had to restrict the size of the initial 
data, which can then been taken arbitrarily large. 

2. In '13t, which studies the case a — b, the author didn't study the Einstein 
evolution equations, which are, as we saw in paragraph 4 the main problem 
to solve. 

3. The present work extends the global existence result established in '24'l, in 
the case a = b and strictly positive cosmological constant A > 0, to the 
case A = 0. 

4-. In the future, we plan to prove the geodesic completness and to relax the 
hypotheses on the collision kernel. 
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